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! Abstract 



We obtain some rigorous results on limit theorems for quantum walks driven by many coins 
introduced by Brun et al. in the long time limit. The results imply that whether the behavior of a 
particle is quantum or classical depends on the three factors: the initial qubit, the number of coins 
M, d = [t/M], where t is time step. Our main theorem shows that we can see a transition from 

Ctf ■ 

^ ' classical behavior to quantum one for a class of three factors. 
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I. INTRODUCTION 



Discrete-time quantum walks are considered as quantum counterparts of discrete-time 
classical random walks . It is often showed that a useful quantum search algorithm 



can be designed based on the quantum walk 



4J. Moreover, the result on a quantum walk on 
Z with an absorbing wall [5j has been applied to solve the transport problems in solid-state 
physics of strongly correlated electron systems in Q|. 

The important factor describing a difference between the quantum walk and the classical 
random walk is the order of variance in the long time limit. Let X t and X t describe the 
position on Z of finding a particle for a symmetric Hadamard walk and a symmetric classical 
random walk, respectively. To get the property appearing to the quantum walk, we will 
consider the limit theorem of weak convergence. The central limit theorem implies X t j\ft =^> 
N(0, 1) as t — > oo, where N(0, 1) is Gaussian with mean and variance 1 and "=>■" means 
weak convergence, while X t /t =>- Z as t — > oo with the following density function (see 

inn 

3, 18|, M for more details): 

p(x) = /( - 1/ ^ 1/ / ^ )(3:) , (1) 

where J^(x) is the indicator function of a set AcK. Thus a quantum particle can travel 
to quadratically farther position than a classical particle on Z. This property also suggests 
the usefulness of quantum walks for spatial quantum searches. 

In this paper, we show a transition from the quantum walk to the classical random walk 
in the view point of the scaling order for the weak convergence with respect to time steps. 
We treat the quantum walk introduced by Brun et al. 10, llll]. If we apply the M coins to 
the quantum walk, the quantum walk is called as an M-coin quantum walk (M-CQW) here. 
The M-CQW is a quantum generalization of the random walk depending on the previous 
M-step memory ^Q]. The quantum coin describing the one step dynamics can be obtained 
by replacing the nonzero entries of an adjacency matrix of the de-Bruijn digraph to some 
nonzero appropriate values. See [3] for more details. Throughout this paper, the initial 
qubit is given by two cases, i.e., (A) the initial qubit described by M-th tensor product of 
if with ||(Ail| 2 = 1 and (B) the initial qubit chosen from the basis of coin space randomly. 
Brun et al. showed the first and second moments of M-CQW in the long time limit 
with the initial qubit (A). We call arcsine law distribution the distribution corresponding to 



2 



arcsine law whose distribution has the following density function: 

s(x) = ItlM. (2) 

7TV 1 — X 2 

We often use the arcsine law distribution to describe the quantum property, for example, 
the distribution appears as the limit density function for continuous-time quantum walks on 
zjjjj, Q and on hypercube It] ]. 

Our main result provides the limit distribution for large time steps given by convolution 
of N(0, 1) and arcsine law distribution on an appropriate setting, and it also provides the 
limit distribution of product of two independent random variables X and Z on another 
appropriate setting, where X has the distribution N(0, 1) and Z has the distribution with 
the density function p(x). 

This paper is organized as follows. Section 2 is devoted to the definition of M-CQW. In 
Section 3, we introduce a useful lemma to obtain the characteristic function of the position 
of finding a particle, and we give a proposition in relation to the weak limit theorems as time 
step t — > oo with t = (mod M) by using the lemma. Section 4 proposes three types to 
find the quantum-to-classical transition and gives the limit theorems. At first, we consider 
the case M < t < 2M for large M. Secondly, we treat the other initial condition, that is, 
a mixture of cases (A) and (B). Finally, we give the weak limit theorems for M = Q(t a ), 
d = e^ 1-0 ) (0 < a < 1) as t -> oo, where we say that 0(/(£)) = g(t) if and only if 
0<hm t ^ oo !/(*)/<?(*) | <oo. 



II. DEFINITION OF M-CQW 

The discrete-time quantum walk on Z is a quantum generalization of the random walk 
with an additional coin state He- The state space is described by direct production of 
Hp (position state) and He (coin state). Here, Hp is associated with standard basis \x), 
x £ Z. He is generated by chiralities corresponding to the directions of the motion of 
the particle. The one step unitary transformation U is described by U = S ■ C with two 
operations, "coin flip" C and "shift" S. The coin flip operates on the coin state by using 
a unitary coin H and leaves the position state. The shift operator moves the particle to 
one unit following the chirality of the particle. Let \l)c — T [l>0]> | — l)c — r [0, 1] , where 
T means the transposed operator. In M-CQW, the standard basis of the coin state is 
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given by \£ M -i)c <8> • • • <8> |£o)c = |Ca/-i, • • • , Co) with £ e {-1, 1}. We prepare M quantum 
coins H , Hi, . . . ,Hm-i, where Hj is given by Hj = j® M -i- 1 g) H ® J®- 7 and if is the 
two-dimensional Hadamard matrix: 



if 



1 

7! 



1 1 
1 -1 



At time t, the coin if mo d[M/] acts on the coin state, where mod[i, M] means remainder of 
t/M: 

C\x, Cm-i, ■ ■ ■ j Cj> • • • j Co) = Cm-i ; • • • > H£j, . . . , £o)- (3) 
where j = mod[t, M\. S moves the particle to one unit following the j-th chirality: 

S\x, Cm-i, • • • , ij, • • • , Co) = k + Cj> Cm-i> • • • ) ij, ■ ■ ■ > Co)- (4) 
We can express U as 

U = S® {f^-j-i Pi + 5-1 g, g, P i g, f ^ 5 (5) 

with 5 = Errezl^ + 1 )( a; l and p i = \j)(j\ H (j = i,- 1 )- Remark that if M = I, the 
quantum walk is equivalent to 2-state Hadamard walk on Z. Assume that the time step is 
described by t = dM + q with < q < M. Let $ [ t M] (x) be the coin state of a particle at 
time t and position x. Let describe the position of finding a particle at time t. The 

probability of finding a particle at time t and position x is defined by 



P(X, 



[M] 



X) 



^\x)\\\ 



The spatial Fourier transformation of &\ M \x) is given by ^ u> (k) = XLez &t U> {x)e lhx . As- 
sume the initial qubit is with ||v?j|| 2 = 1- Then from Eq. ([5]), we have 



J M ]/ 



M 



W\k) 



<g) *T\k)®<$^l{k), 

j=g+l j=l 



(6) 



where V { p\k) = H d (k)y j with H(k) = (e ik \l)(l\ + e~ ik \ - 1)(-1\)H. 
Miyazaki et al. 19| described M-CQW by their quantum walk with a quantum coin de- 
scribed by the (2 7 + l)-dimensional unitary representation of the rotation operator with 
half-integer j. The M-th tensor-product of C 2 space Wm is decomposed into (2 7 + 1)- 
dimensional irreducible rotation group modules Vj with half-integer j; Wm — (Bjdf^Vj 
with 

M \ ( M \ 
(M-2j)/2j ~ V(M-2j)/2-lJ- 



d [M] 



The multiplicity dj is obtained by using the highest weight decomposition (see p. 66 in 
[20I ) . Thus the limit density function of M-CQW is described by a linear combination 
of their density functions. However it seems to be complicated to compute the unitary 
basis transformation operator Km '■ VVju — * d^Vj, so in this paper, we treat M-CQW 
directly. 

III. THE CHARACTERISTIC FUNCTION FOR M-CQW 

We consider the following two cases for the initial qubit $0; 

Case (A): % = ^ M , with y? = T [l/V2,i/V2}. 



Case (B): $ = <f>M-i ® ' ' ' ® 0o? where {0j} J=o is an i.i.d. sequence of Bernoulli 
random variables with Pr(0j = 1) = Pr (</>•,■ = —1) = 1/2. 



Remark that the initial qubit of Case (A 



is in a pure state, while the initial qubit of Case 
(B) is in a mixed state [21[ . Brun et al. H[ computed the first and second moments of M- 
CQW with the initial qubit of Case (A) in the long time limit and showed that the variance 
grows in proportion to square time steps. By definition of the spatial Fourier transformation 
of ^> d M \x), we have 

E(e«xr) =jT<*r(*), 4r<*+ o>|- pi 

Then applying Eq. ([6]) to the above equation, we can provide the following lemma with 
respect to the characteristic function for x\ M ^ to obtain some weak convergence theorems: 

Lemma 1. Let Q d (k,£) = (¥f \k),¥f o) (k + ^)) and C d {k, f) = Tr[H- d (k) ■ H d {k + £)]/2. 
Assume that t = dM + q with < q < M. Then the characteristic function of x\ M ^ can be 
expressed as 



Case (A): 



Case (B): 



E (e^) = J" Q d+1 (k,O q Q d (k,0 M ~ q ^- (8) 
E (e^ MI ) = J^C d+1 (k,O q C d (k,0 M -^. (9) 
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By Eq. ©, we have $[ M] (x) = J2 VM +- +m = x P vm ® " " ' ® for t = M. Note that 

(Pii],PjT]) = Sij/2 for r] G {00)1 — 1); Therefore xj M ' has a binomial distribution 

B(M, 1/2). 

Proposition 1. For both Cases (A) and (B), X t /y/i =^ N(0, 1) as t — > oo with t < M. 

In the rest of this section, we will show a weak limit theorem taking t — > oo with g = 0, 
i.e., t = dM. Put 

which is the density function for the symmetric Hadamard walk [7|, la, |9j in the long time 
limit. 

Proposition 2. Assume that t = dM . 

1. If we fix M and take d — > oo, then with the initial qubit of Case (A) (resp. (B)), 
x\ M] /t (resp. T7te limit density functions /M fresp. ^Mj o/yW 

(resp. Z' M 1 j zs gzi>en fry; 
Case fAj: 

[M/2] /M\ ofr 1 
= J2 { j ) \l- 2 j/M\ V{Xj)M ~ jQ{Xj)j + c V=-n) (10) 

with xj = x/\l - 2j/M\, 

V{x) = (l + Vl - 2x 2 ) /2, Q(x) = (l - Vl - 2x 2 ) /2, (11) 

where c is determined by 



oo 

f [M \x)dx = 1. 



Case 

A/ [M/2 



9 [M| M = (i) £ C') jl^M, + <V-«.,4(<fa), (12) 



i=0 

where Xj = x/\l — 2j/M\ and d is determined by 

gi M \ X )dx=l. 
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2. If we fix d > 2 and take M — ► oo ; £/ien X^/t Yd with the initial qubit of Case (A), 
while Xl M] /Vi Z d with the initial qubit of Case (B). Let p d (k) = (¥ d vo \ D k ¥f o) )/d, 
and i>d{k) = Ti[H~ d (k) ■ D\H d {k)]/2d with D k = id/dk. The limit distributions Y d 
and Z d are given by 
Case (A): 

P(Y d <x)= [ * (13) 



Case (B): 



f f e~ u2 l 2 dudk , . 

P(Z d <x)= / . (14) 

J ^{(fc,«)6[0,27r)xffi:«(^(fc)) 1 /2< 2: } [Alt) ' 

To prove Proposition 2, we will use the following lemma. Let the eigenvalues and eigenvectors 
of H(k) be denoted by Xo(k), Xi(k) and \vo(k)), \vx(k)), respectively. 

Lemma 2. 1. Let hj(k) = D k (Xj(k))/Xj(k), (j G {0, 1};. Then 

hoW + h^h) = 0. 

2. Define p(k) = \(v {k),(f )\ 2 and q(k) = 1 - p(k) = \(vi(k), if )\ 2 . Then 

p{k)+q(k) = 1. 



3. Let k(x) = arccos(x/v / l — x 2 ). Then 



2?r jt, poo 







dk 

g(hj(k))w(k)— = j g(x)w(k(x))p(x)dx. (jG {0,1}) 



From now on, we will prove Proposition [2j 



Proof of Proposition [2l 

1. Let Xj(k) = e ie ^ (j = 0, 1). We should remark that (0j(k + £/d) - Bj{k))d = 
£hj(k) + Oid- 1 ) and \ Vj (k + £/d))(vj(k + = \vj(k))(vj(k)\ +0(d~ 1 ). By Lemma 
E]© and P), we have for fixed M, 

Q d {k,£/t) = e^M^+or 1 ))/^^) + g-i^M^or 1 ))^^) + o^-i^ (15) 
C d (jfc,£/t) = i ( e *(M*)+o(t-i))/Jtf + e-^MfcJ+oct- 1 ))/^^ + 0(r i ); (16) 
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where h(k) = ho(k). Combining Lemma [T] with Eq. (fT5|) . the characteristic function 
with the initial qubit of Case (A) is described as 

;/ / w\ r 2 * . ... _ ,//, 

(17) 



i=o v J 7 



27r e ^(i- 2 ,/Af)Mfcyf- J(A;)g , (A , ) g_ 



From Lemma [2] ([3]), we obtain the desired conclusion. We can also prove Case (B) by 
combining Lemma [1] with Eq. (fT6j) . 

2. Noting that £>£#(£;) = a 3 H(k), we have Tr(H~ d (k) ■ D k H d (k)) = for all d, where cr 3 
is the Pauli matrix: 

r i 



0-3 



Then we obtain for fixed d, 



-1 



C d (k,at e ) = l-^ d (k) + 0(t-n, 
with 6 1 > 0. By applying Eq. ( IT91 to Lemma [H, in Cases (A) and (B), we see 



(18) 
(19) 



lim E 
lim E 



Ax\ M] /t 



2- 



dk 



2tt : 



2tt 



2tt' 



(20) 
(21) 



respectively. Thus Eq. ( 1201) gives the desired conclusion in Case (A). For Case (B), 
noting that e~^ 2 ^ 2 is the characteristic function of N(0, 1), Eq. f|2Tj) can be rewritten 

as 



lim E 


~ e izxM/vi 


p2ir r 






Jo J- 



dudk 



□ 



In the case of the initial qubit \\)® M for Proposition [2] (1) and (2), we can obtain 
£[(y[M])2j = x _ 5 /( 4v /2) + l/(4M v / 2) and E[Z^_= 1/8 (d = 2), 7/72 (d = 3) in a 



similar way. The results were shown in Brun et al. [11]. 

In Proposition [2] (2) (d = 2 case), we have /^(A:) = (sin2/c)/2, /^(^O = 1- Therefore we 
obtain the following corollary. 
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Corollary 1. As t -»• oo with t = 2M, in Case (A) (resp. (B)), X^/t Y 2 (resp. Z 2 ), 
where Y 2 has a scaled arcsine law distribution (resp. N(0, 1)). The limit density function of 
Y 2 is given by 

2s(2x) 



2/(_i /2 ,l/2)(z) 



tx\/\ - 4x 2 
where s(x) is defined by Eq. (TJ|). 



IV. A CROSSOVER FROM CLASSICAL BEHAVIOR TO QUANTUM ONE 

As Brun et al. claimed in 1101, the behavior of M-CQW with the initial qubit of Case 

nn 

(A) remains "quantum" in contrast with [17|, ll8j in the long time limit. However by 
considering the initial qubit of Case (B) and time steps t ^ (mod M), we can see a 
"classical" property. From Propositions Q] and [2] (j2J) with the initial qubit of Case (A), the 
behavior of a particle becomes classical as t — > oo with t < M, while the behavior grows 
quantum with t > 2M in the view point scaling order. So we will consider long time limit 
with M < t < 2M. By Proposition [2] ([2]), the scaling order of the weak convergence is given 
by t and \/t in Cases (A) and (B), respectively. To find a behavior corresponding to a kind 
of quantum to classical transition, we will introduce another initial qubit, that is, a mixture 
of Cases (A) and (B). From Proposition [2] (TjQ) and (j2J) with the initial qubit of Case (B), the 
scaling order grows t for fixed M, while it grows y/i for fixed d in the long time limit. So 
we analyze the limit theorem for d, M — > oo simultaneously with the initial qubit of Case 

(B) . Therefore, we consider the following three assumptions. Let < (3 < 1. 

Assumption (a), t = M + M 13 with the initial qubit of Case (A). 

Assumption (b). t = 2M with a mixture of initial qubits of Cases (A) and (B), i.e., 
$ = Lp® M>3 (g) 0i (g> • • ■ <8> 4> M _ M p, where {4>j}jLl Ml3 is an i.i.d. sequence with <j>j = |1), |— 1) 
with probability 1/2, respectively. 

Assumption (c). M ~ t 1 ' 13 , d ~ t 13 with the initial qubit of Case (B), where f(x) ~ g(x) 
means lim^oo f{x)/g(x) = 1. 

We obtain a phase diagram in relation to the limit distribution in each case; Define 
D = {((3,6) G [0,1] 2 : 9 > max{l/2, /?}}, D' = {((3,6) G [0, l] 2 : 6 > (l+(3)/2}. Let the 
scaled arcsine law distribution with the density function 2 /3 s(2 f3 x) be denoted by 
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Theorem 1. Let X, Z, and be independent random variables with X ~ N(0, 1), 

Z ~ K , and ~ where the distribution K has the density function p(x) defined by 
Eq. (QP. As t -> oo ; 

1. Under Assumption (a), 



x\ M] /t b 



S (x) 


if (13,6) eD, 


X 


if 9 = 1/2 and < < 1/2, 


X + W (0) 


if 9 = (3 =1/2, 


w (0) 


if 1/2 < 9 = (3 < 1, 




if 9 = (3 = 1. 



2. Under Assumption (b), 



X\ M] /t 6 



5 (x) 
X 



if((3,9)eD, 

if 9 = 1/2 and0<(3< 1/2, 



3. Under Assumption (c), 



x\ M] /t 



X + W^ if9 = (3 =1/2, 
W (/3) i/ 1/2 < 9 = (3 < 1. 

6 Q (x) if{(3,9)eD', 

X if 9 = 1/2 and (3 = 0, 

XZ ifO < 9 < 1/2 and 9 = (1 +(3)/2, 

Z if 9 = 1 and 13 =1. 



Under Assumptions (a) and (b), as long as f3 < 1/2, the behavior is classical, while if (3 > 1/2, 
then the standard deviation of a particle is proportional to t@ and its limit distribution is 
described as a scaled arcsine law distribution. The cross point of the two regions, that is, 
(3 = 9 = 1/2, the limit distribution is obtained by the convolution of N(0, 1) and an arcsine 
law distribution. Under Assumption (c), the classical (resp. quantum) behavior appears at 
the end point (3 = (resp. (3 = 1). If < (3 < 1, the scaling order of the weak convergence 
is t( 1+ ^)/ 2 and its limit distribution is expressed as product of two independent variables, X 
and Z. 
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Proof. 1. Lemma [I] implies that the characteristic function for X t /t e under Assumption 



(a J is 

r2vr 



B ( e rfV) = jf ^(wO^^AT^- (22) 

Remark that m(k) = 0, Qi(k,£/t e ) = 1 - £ 2 /( 2 + 0(t~ 3e ). From the estimation 



and Eq. (fT8|) . as t — > oo, 



1 if ^ > 1/2, 

-< (23) 
e"« 2 /2 if e = i/2. 



.1 if 9 < (3, 

Q 2 (k,Z/tY^{ (24) 



When 9 = (5 = 1/2, we see that by applying Eqs. and ([24]) to Eq. (J22 



oo «/ — oo 



V 27T 



since e~^ 2//2 is the characteristic function of N(0, 1) and /^(&) = sin2/c. For other 
cases except 9 = (3 = 1, the desired conclusion can be derived similarly. Moreover 
Corollary CD (Case A) gives the result of 9 = (3 = 1. Therefore, we complete the proof 
of part (1). 

2. The characteristic function under Assumption (b) can be expressed as 

E ( e ^rv^ = J^ 2 (k,atT c,(k,atT- M0 ^- (25) 



By Eq. QT9J), we get 

|l if > 1/2, 

C a (*,£/0-< (26) 
[e~« 2 / 2 if 0=1/2. 

From Eqs. ( l25l) and ( l26l) . we have the desired conclusion in a similar fashion of part 
(!)• 

3. Under Assumption (c), Lemma [1] (Case (B)) implies 

E (f$«l*)= j^C d {k^r d ± (27) 
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Eq. ([16]) yields 

C d (k^/t e ) = 1 - ^h 2 (k)t 2 ^+o(t^-% 
in the condition of (3 — 9 < 0. So we have as t — > oo, 



1-/3 



,1 if 0> (l + /3)/2, 

{(^(A.e/O} H (28) 

e -W*)/2 if0 = ( 1 + /S )/ 2 _ 



Note that 



e i€M*)«i^ dw . (29) 

-oo V27T 

Applying Lemma [2] © and Eqs. (T58D d29j) to Eq. (T57D, we have for = (1 + /3)/2, 

-ooJ-oo V27T 

Proposition [T] and Proposition [2] (Op) for M = 1 give results of (/3, 0) = (0,1/2) and 
((3, 9) = (1, 1). Therefore we complete the proof of part (3). 

□ 



Under assumption t = dM, Eqs. ( 1161) and ( 1191) imply that for < (3 < 1 
lim{c>(MA-)} 

= lim lim <Cd\k, 



1-/3 r / p \ ^ M 

lim lim | fc, 



Af->oo d^oo [ y yMd 

oo Af^oo [ " V ' VMrf 
£ 2 /i 2 (fc)/2 



e 



Similarly, from Eqs. ( fl5l) and ( 1T81) . we see that for < (3 < 1, 



lim {Q t0 (k,£/t)} = lim lim <Q d [k, 

a. 



M 



M 



= lim lim < Qd [ k, , , , 

d-ooM^oo y ' Md 

_ e «fc(*)(P(*)-9(*)) . 

Therefore, applying the above equations to Lemma [H we get the following result. 
Proposition 3. Let £ = c?M. VKe impose the same assumtion in Proposition^ 
1. y/MZW XZ (M ^oo), and Z d =>- XZ (d oo). 
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2. Assume that d ~ t 13 and M ~ t 1 13 with initial qubit of Case (A). Then we have 
A"j M ' jt =>- Y as t — > oo, for < (3 < 1, where Y has the density function 

fix) = 3I CV^/Vs)( x "> 
tt(1 + x 2 )V1-8x 2 ' 

Furthermore, 7' M ' =>- Y (M ^ oo), and Y^^Y (d—>oo). 

In the case of the initial qubit |1}® M for part (2), we can obtain the limit distribution Y 
similarly, where Y has the density function 



fix) 



7r(l - X)yj(l ~2x)x' 



This gives £"[y 2 ] = 1 — 5/(4y / 2) which was shown by Brun et al. 11] 
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